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Abstract 

The possible paralelism existing between phase transitions and fracture in disor- 
dered materials, is discussed using the well-known Fiber Bundle Models and a 
probabilistic approach suited to smooth fluctuations near the critical point. Two 
limiting cases of load redistribution are analyzed: the global transfer scheme, and 
the local transfer rule. The models are then studied and contrasted by defining the 
branching ratio as an order parameter indicative of the distance of the system to the 
critical point. In the case of long range interactions, i.e., the global rule, the results 
indicate that fracture can be seen as a second-order phase transition, whereas for 
the case of short range interactions (the local transfer rule) the bundle fails suddenly 
with no prior significant precursory activity signaling the imminent collapse of the 
system, this case being a first-order like phase transition. 
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1 Introduction 



Fracture phenomena have attracted a lot of interest in the last several years, 
both experimentally and theoretically. In the lab, a disordered material sub- 
jected to an increasing external load can be experimentally studied by mea- 
suring the acoustic emissions before the global rupture. It has been shown 
[1,2] that this intense precursory activity in the form of bursts of different 
microscopic sizes follows a well-defined power law. From the theory side, the 
understanding of fracture in heterogeneous materials has progressed due to 
the use of lattice models and large scale numerical simulations [3] . Recently, 
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the introduction of models of material failure has lead to the evidence that 
rupture can be viewed as a kind of critical phenomenon [3-6]. Nevertheless, 
the question of whether rupture exhibits the properties of a first-order or a 
second-order phase transition remains under discussion as well as what is the 
order parameter that determines the type of transition. 

In this field, it is important to use models able to describe the complexity of the 
rupture process, although they should be simple enough to permit analytical 
insights. To this class of models belong the well-known Fiber Bundle Models 
(FBM) widely used since their introduction more than forty years ago [7,8]. In 
static FBM, a set of fibers (elements) is located on a supporting lattice and one 
assigns to its elements a random strength threshold sampled from a probability 
distribution. The lattice is loaded and fibers break if their loads exceed their 
threshold values. Now, one can assume different load transfer rules to mimic 
the range of interactions among the fibers in the set. The Global Load Sharing 
rule is the simplest theoretical approach one can adopt to make the problem 
analytically tractable, which implies that the load carried by failed elements is 
equally distributed among the surviving elements of the system, representing 
in this way a long-range interaction among the constituent parts of the system. 
This is a kind of mean field approach to the more extreme Local Load Sharing 
scheme, where the interaction among the elements is short-range and the load 
supported by the failing elements is transferred only to a neighborhood. 

Very recently [9] , we have developed a probabilistic approach suited to smooth 
the fiuctuations around the point of final collapse. In this paper, we extend 
the method to analyze the local load sharing scheme and contrast the two 
limiting cases of range interactions: the long range and the short range. While 
for the case of global load sharing, we obtain different scaling relations that 
point out that the fracture of a fiber set with long range interactions may be 
a continuous transition, for the case of short range interactions, the failure 
exhibits no significant precursory activity and the system breaks abruptly 
showing, in this case, the properties of a first order transition where no scaling 
laws can be identified. The rest of the paper is organized as follow. In Section 
2, we present the standard approach to the FBM. In Section 3, we introduce 
the new probabilistic method and the results derived from it and from Monte 
Carlo simulations. Finally, Section 4 is devoted to discussion and conclusions. 



2 Static Fiber Bundle Models 



Let us first recall the basic ingredients of the static FBM and how one proceeds 
in numerical simulations. The term static means that time plays no role in the 
model [10]. The system under consideration is a set of Nq elements (fibers) 
located in a supporting lattice each one having at the initial state a zero load 
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and a fixed strength threshold vahie (quenched disorder) sampled randomly 
from a probability distribution P{cr). The system is then subjected to an 
external force F which is shared democratically among the fibers, that is, the 
distribution of the total force produces that each element increases its load, 
(7, in the same amount. This individual stress acts as the control parameter. 
The loading process is done quasistatically; i.e., the external force is increased 
at a sufficiently slow rate as to produce a single breaking event when the 
stress in the weakest element equals its threshold value. Then, the increase of 
external force stops and the load of the broken fiber is transferred according 
to the transfer scheme assumed. In the global load sharing FBM, this impfies 
that the load on fiber i is given by cTj = F/ns{F), where Ug is the number of 
surviving elements for a given external load F. In the local load sharing FBM, 
the above relation does not hold any more because there the redistribution of 
load is performed among the nearest neighbors of the failed fiber giving rise 
to the appearance of regions with different stress concentrations. 

In both schemes, the rupture of a ffber may induce secondary failures which 
in turn may trigger more failures. This process of induced failure at constant 
external load, termed avalanche, stops when all surviving elements carry a 
load lower than their thresholds. The system is then loaded again and the 
process is repeated until the ffnal catastrophic avalanche provokes the total 
rupture of the material, which occurs at a critical load o"c that depends on 
the probability distribution from where the individual strengths were drawn 
as well as on the system size. At this point, it is worth recalling that the 
exact value of Uc can be analytically obtained in the thermodynamic limit 
for the global load sharing case while for local load sharing schemes there is 
no theoretical approach leading to dc. It is also known that the critical load 
(Jc for the global load sharing FMB is non-zero and independent of A^o for 
Nq ^ oo, whereas cXc tends to zero as Nq goes to inffnity in the case of local 
load sharing models [8,11]. Fiber bundle models have also been recently used 
in self-organized criticality (SOC), a theoretical framework widely used for the 
study of avalanche phenomena in disordered systems. It has been shown using 
these models that systems with plastic behavior can reach a SOC state just 
before the global rupture [12]. A second case of self-organization with power 
law distributions in several quantities corresponds to the situation in which 
the fracture process coexists with a healing process [13]. 

In numerical simulations, the cycle of complete breakdown of the material 
is performed many times in order to average out the effect of ffuctuations. 
Nevertheless, the stress history of a particular element is made of steps, so that 
the fluctuations around the critical load can never be completely avoided. Let 
fiber k be supporting a stress o"^ at a given step of the failure process. It will 
continue to support until a global driving or a stress transfer from one (or 
more) of the failed elements occurs. At this moment, element k instantaneously 
changes its load to al — + adist, where adist is the load coming from the 
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failed fibers or from the global driving. So, in a later step, clement k could 
receive again load, suffering a second step-like increase in stress. This step-like 
stress history continues until fiber k fails. As we are interested in studying 
the behavior of the system as the critical point is approached, it is of utmost 
importance to find a simple method able to capture the evolution of the system 
near the critical point avoiding as much as possible the fiuctuations. 



3 Probabilistic Approach to Static Fiber Bundle Models 



Let there be a set of Nq elements located on a supporting lattice. Suppose that 
each element carries a given load a, which is set to zero at the initial state. 
Fibers break depending on their strengths which are distributed according 
to a probability distribution P{a). Different probability distributions can be 
considered. In materials science the WeibuU distribution is widely used, 

P(a) = l-e"^^^', (1) 



where p is the so-called Weibull index, which controls the degree of disorder 
in the system (the bigger the Weibull index, the smaller the disorder), and cto 
is a load of reference. In the following, we will assume ctq = 1, and therefore 
the loads are dimensionless. Although the results we show hereafter are for 
the Weibull distribution to gain in definiteness, they have been also obtained 
for a wide class of distributions. Besides, for continuous distributions decaying 
fast enough, the approach to the critical state does not depend on the detailed 
form of the disorder [14] . 

Equation (1) represents the probability that an element fails under the indi- 
vidual load a. Now, consider the case in which an element drawn from Eq. (1) 
supports a load cxi, but breaks under a new load a2- The probability that this 
happens is given by 

, P(CT2) - P(ai) 



which is equal to 



p((7i,(72) = l-e-(-^-^). (3) 



So, the probability q{ai, (T2) that an clement that has survived to the load ai 
also survives to the load (72 will be given by 

, (72) = 1 - , £72) = e-^'^^ -<) . (4) 
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All these probabilities depend on the state of stress of each element, which is a 
complex function of both the control parameter and the stress redistributions 
due to fiber failures. To mimic the quasistatic increase in load on the system as 
is applied in MC simulations, we impose the condition that under an external 
force F, the next breaking event consists of one single failure. Now, we will 
explore how the driving process controls the approach to global failure in the 
global load sharing case and in a particular version of the local load sharing 
schemes. 



3.1 Global Load Sharing case 



Let suppose that after the latest avalanche, there are Nk surviving elements 
each one bearing a load 0"^. The assumption of quasistatic increase in load on 
the system implies that the load ai needed to provoke the failure of only one 
element is [91, 



(7i 



(5) 



where in Eq. (5) Nk = No and cTfc = at the initial state. Elevating the external 
force up to the N^ ■ (Ti level a first element breaks. As we are dealing with a 
global load sharing set, the choice of the broken element is irrelevant because 
all of them are equivalent. Once the first element fails, the redistribution of 
its stress takes place which may induce secondary failures and so on, until the 
end of the avalanche. 

Now, the number of surviving elements in the set after an avalanche has come 
to an end can be calculated by the recursive relations [9], 

Nj+i^ Nj-q{aj_i,aj), (6) 

Nj ■ C7j = Nj^, ■ (7,_i (7) 

Nj = Nj+, (8) 



Equation (6) relates the number of surviving elements between two successive 
avalanche steps, whereas Eq. (7) allows the computation of the new load acting 
on the Nj surviving elements. This is possible because we are under a global 
load sharing scheme and thus all surviving elements support the same load. 
Finally, Eq. (8) determines the end of the ongoing avalanche. 

The dynamics of the system is completely determined by Eq. (5), (6), (7). In 
this way, the size of an avalanche is given by the number of elements that break 
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between two successive steps of external loading. The total stress accumulated 
in the system can be calculated multiplying the number of intact fibers before 
an avalanche starts by the load given by Eq. (5). The critical load, defined 
as the load needed to provoke the total collapse of the system, is equal to 
the load on the intact fibers just before the final catastrophic avalanche. Note 
that in this probabilistic approach, in contrast to Monte Carlo simulations, we 
need to store only the information concerning the loads of the intact elements, 
that is, the threshold dynamics is omitted with the subsequent advantages of 
saving computer resources and the possibility of exploring systems of larger 
size. 

For Eq. (8) to be satisfied, we can proceed in two different ways in order to 
determine when an avalanche ends, to which we will refer as the continuous 
and the discrete cases. For the continuous case, the number Nj^i of surviving 
elements is considered a real number. This means that condition (8) is never 
fulfilled before the final avalanche. So, condition (8) is replaced in numerical 
simulations by a factor z/ <C 1 that determines the end of an avalanche, i.e., 
if Nj — iVj_|_i < u the avalanche stops; otherwise it continues. In the discrete 
case, Nj^i is considered to be an integer number, so that after each iteration 
of Eq. (6), Nj^i has to be rounded up. This is done comparing the remainder 
of Nj^i, A, with a random number a uniformly distributed in the interval 
[0, 1[. Thus, if q; > A, A^j+i is equal to its whole part, and if not, A^j+i is 
equal to its whole part plus one . Next, we check whether condition (8) is 
satisfied for the rounded value of Nj^i or if a new iteration of Eq. (6) has 
to be performed. The continuous approach has the great advantage that the 
fluctuations arc ruled out, whereas for the discrete case the results are similar 
to those obtained by Monte Carlo simulations where it is necessary to average 
over many realizations in order to get accurate mean values. This is because 
in the global fiber bundle model the central limit theorem applies [15]. 

Now, we proceed to explore the behavior of the system near the critical point. 
In particular, we are interested in inspecting the evolution of some quantities 
as the critical point is reached. In order to avoid unnecessary fiuctuations 
around the critical point, the results shown below have been obtained for 
the continuous case of the probabilistic approach (p = 2). It turns out [9] 
that the avalanche size near to the critical point diverges with an exponent 
7 = I as s ~ ((Tc — <^)~'^- A similar behavior, through a mapping of a fuse 
network model to the global fiber bundle model used here, has been recently 
reported [5]. This mapping between fiber bundle models and fuse networks 
with strong disorder was first noted [16] a few years ago; nevertheless, it has 
also been pointed out that 3D fuse networks apparently do not follow the FBM 
picture [17]. Note, additionally, that for the probabilistic continuous version, 
this relation is obtained for a single realization avoiding, in this way, the large 
number of iterations performed in MC simulations. In Fig. 1 we have depicted 
the derivative of the number of broken fibers, (dN/da), as a function of the 
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distance to the critical point cTc — cr in a log- log plot. This rate dN/da diverges 
as ((Tc — (t)"'*' also with 7 = |, thus qualifying a critical mean field behavior as 
was already shown [6] by means of analytical analysis of fiber bundle models. 
In Ref. [12], a similar scahng behavior is addressed for the derivative of the 
strain carried by the fibers with respect to the driving force. 

Another way to shed light on the critical behavior of this type of system is to 
define a branching ratio C, for each avalanche. This magnitude represents the 
probability to trigger future breaking events given an initial individual failure 
[18,19], and is related to the number of broken fibers by 

C=ll^. (9) 

< 2; > 



The above relation can be obtained by thinking of the evolution of fracture 
as a kind of branching process [20]. In this process, each node gives rise to a 
number n of new branches in the next time step. The average number < n > 
of new branches is called the branching ratio. Let us denote by nt the number 
of branches at a given step t of the branching process, and by t^ax the total 
number of time steps before it dies. Then, 



t-m ax 1 

E rit+i 
t=o 

+ : 

E nt 
t=o 



and 



no 



t rna x 

E n-t 
t=o 



As no = 1, C = 1 ~ ~7 where ntot is the total number of nodes developed 
in the branching process. For a fracture process, ntot is equal to the average 
number of failure events. So, Eq. (9) defines the branching ratio. We represent 
by < 2; > the average number of elements that fail in one avalanche, which 
is a function of the control parameter a and coincides with s. This analogy 
between fracture and branching processes has been previously used to study 
the criticality in the process of fragmentation of Hg drops [21]. The branching 
ratio will then act as the order parameter. It takes the value 1 when the system 
is critical thereby representing a measure of the distance of the system from 
the critical state [19]. We would like to remark here two characteristics of the 
branching ratio defined as above. First, it coincides with its general definition 
for a branching process. Here, the values of the branching ratio are always 
less than or equal to one because we are analyzing an irreversible breaking 
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process that cannot continue forever, so that a value of the branching ratio 
greater than one would imply a physically unreachable situation. Secondly, we 
have introduced a definition in terms of the avalanche size and not in terms 
of the cluster size. The avalanche size in fiber bundle models is a measure of 
causally connected broken sites while the cluster size is a measure of spatially 
connected broken sites. It is clear that in the global fiber bundle model, the 
spatial correlations are ruled out and then the distribution of avalanche sizes 
does not coincide with the distribution of cluster sizes. 

We have numerically computed C for the continuous version of the probabilistic 
model. The results obtained for a system of A'^o = 50000 elements and p = 2 
have been plotted in Fig. 2. It can be seen in this figure that the branching ratio 
approaches the unity as the critical load is reached. In the figure, the values 
of C are collected for all the avalanches except for that which provokes the 
catastrophic event leading to the collapse of the system. So, this dependency 
of C with a means that very close to the critical point, the initial failure 
provokes the breaking of another fiber which in turn induces tertiary ruptures 
and so on. This is the result that should be expected from the divergence 
of the avalanche size at the critical point. Figure 3 shows the behavior of the 
order parameter near the critical point. As can be seen, near the critical point, 
the relation 1 — ^ ~ (o^c -~ o")'^, where /9 = | applies. Note the similarity of 
the behavior with those obtained for the magnetization in known magnetic 
systems with second-order phase transitions. It is of additional interest to 
note that the branching ratio also captures the feature that the precursory 
activity is significant only for strong disorder (for the WeibuU distribution, 
small values of p) ; whereas for small disorder the system behaves more similar 
to the breakdown of homogeneous materials. On the other hand, the branching 
ratio does not depend on the size of the system for large systems, in contrast 
to previous results in other fracturing systems [19]. In Fig. 4 we illustrate 
this behavior by plotting the value of the branching ratio for the last non- 
catastrophic avalanche, C*, as a function of the system size. In all cases, the 
numerical simulations were performed for p = 2. It is clear from the figure 
that (* reaches the unity as A'o goes to infinity. 



3.2 Local Load Sharing Case 

It is well-known that the local load sharing FBM is much more complicated 
than the global load sharing case. The complexity of the fracture problem 
increases because the load borne by failed elements is transferred to nearest 
neighbors and then there appear regions of stress concentration throughout 
the system. The distribution of load is now not homogeneous, and we have to 
carry on several lists to record the individual load of the fibers in the system. 
However, with the probabilistic approach, we are able to study a few things 
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about the behavior of the system, which are sufficient to remark the great 
differences between the long range and the short range interactions schemes. 
One important difference between the two schemes is that contrary to the 
global load sharing case, in the local scheme there is no significant precursory 
activity signaling the approach to the final collapse and the system undergoes 
an abrupt catastrophic avalanche. For maximum simplicity, we consider next 
a one-dimensional periodic local load system where the stress transfer is done 
by adding the load of the failed element always to the element on its right. 
Despite of its simpUcity, it has been shown [22] that the general properties of 
this particular model are identical to those of more complex local schemes. 

The probability that the system fails in only one avalanche is given by: 

p((7i, (72) • p((7i, (73) ■ p((Ti, 0-4) . . .p((7i, (7iVo), (10) 

where (jj = iai and ai is 

That is, the probability of having a one-step failure is given by the probability 
that an element fails under a load (T2 having survived to the load cti, multiplied 
by the probability that the next element in the lattice also fails under the new 
load as having survived to the load ai and so on. 

On the other hand, Eq. (4) takes the form 

g(i,j) = e-<(^"-''), (11) 

which gives the probability that an element that has survived to the load i ■ <Ti 
also survives to a new load j • ci, being ci the load that produces the first 
breaking event in the system at the initial state. The probability that the first 
avalanche consists of, for example in one failure, is given by g(l, 2). The first 
avalanche will have size two with probability p(l,2) ■ g(l,3). In general, the 
probability that the first avalanche does not provoke the rupture of the whole 
system, i.e., be finite, is given by 

00 

Pi = ^n;-g(i,/ + i), (12) 
1=1 

where 

I 

Ui = p(l, 2) • p(l, 3) . . .p(l, = l[p{l, i), l< No, 

i=2 
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with III = 1. Thus, the average size of the first avalanche, derived from Eq. 
(12) is 

oo 

S, = ^l-Ui-q{l,l + l). (13) 
1=1 

Numerical simulations of the continuous case, the discrete approximation 
(both defined as above for the global load sharing scheme) and Monte Carlo 
method confirm the validity of Eq. (13). For example, for p = 2 and Nq = 1000 
we obtain 5*1 = 1.0030, 5*1 = 1.0022 and = 1.0039 respectively. This result 
is valid only for the first avalanche in the system. 

We have also verified by means of Monte Carlo simulations that for this lo- 
cal transfer rule, there is no power law in the distribution of avalanche sizes. 
Besides, no scaling relations appear and the system is more sensitive to pa- 
rameters such as the degree of disorder (p) and the size of the system A^o- 
Catastrophic avalanches arise very often in the first stages of the rupture pro- 
cess. 

Finally, we illustrate how one can describe the precursory activity when it is 
limited to a few steps before the global rupture. Let us solve the particular 
situation in which the first avalanche is finite and the second provokes the 
final breakdown of the system. The load needed to break the first element at 
the initial state is 

g(l,2) ^e"^''-'^'^. (14) 

If the value of g(l,2) in Eq. (14) is bigger than 0.5, then the fracture will 
hkely consist of one single event of rupture and the avalanche will not progress 
beyond this first failure. As shown above, for p = 2 a system size of Nq = 1000 
elements is enough to produce first avalanches of size 1 on average. Next, the 
driving force is increased and the new load acting on the surviving A'o — 1 
fibers is 

that is, = In the derivation of o"2 it has been assumed that the second 
avalanche starts with high probability in a fiber different of that which is 
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located just in the crack tip of the first avalanche. This approximation is well- 
justified by numerical simulations of this local load model for which it is known 
that the first breaking events are randomly distributed in the sample, leading 
to the appearance of several small cracks which much later coalesce and grow 
provoking the final breakdown. Now, there are (A^o — 2) elements bearing a 
load (72 and one element supporting a load o\. The system will experience 
on average the final breakdown in this second breaking event if the probability 

g'(l,2) = e-^"'-')-^ (15) 

is lower than 0.5. Nevertheless, ^(1, 2) and ^'(1, 2) are not independent since 
the relation 

g'(l,2) = (g(l,2)f (16) 

has to be verified. This dependence, together with the restrictions g(l, 2) > 0.5 
and g'(l,2) < 0.5 leads to 0.5 < g(l,2) < 0.7 and 0.25 < g'(l>2) < 0.5 as 
acceptable values for ?(1,2) and ^'(1,2) in order to have a two-step system 
collapse. This approximation seems reasonable taking into account the way in 
which numerical simulations of the continuous and the discrete versions of the 
probabilistic approach proceed. So, the condition for the first avalanche to be 
finite and that the second breaks the system is reduced to 

A^o ^ 2''+^ - 2 (17) 

for (7(1, 2) = 0.6, a value in the middle of its range. This relationship between 
A^o and p is confirmed in numerical simulations. The situation considered above 
corresponds to a very brittle rupture since the ratio of the size of the typical 
local damage that induces the system failure to the system size is very close to 
one (iVo — l)/A'"o. Expression (17), although very simple, gives a rough estimate 
of the qualitative relation between the size of the system and the amount of 
disorder when the degree of brittleness of the fracture is high. The situation 
analyzed becomes important and more accurate as we go to larger system sizes 
because for the local load transfer rule and a larger number of fibers, the 
probability to find a weak region somewhere in the system is high and due to 
the local redistribution rule the crack wiU propagate fast enough to break the 
whole system in a few steps (strictly speaking, in the thermodynamic limit 
the critical load is zero) . 
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4 Discussion and Conclusions 



We have proposed a probabilistic approach to fiber bundle models of fracture. 
The cases of long range interactions among the fibers of the bundle and a local 
load sharing transfer scheme (short range interactions) were considered in or- 
der to investigate the behavior of the system near the global breakdown of the 
material. The results obtained for the local scheme indicate that the system 
undergoes a kind of first-order phase transition in agreement with previous 
reports [5,23]. The system fails with no significant precursors announcing the 
incipient rupture, and no scaling relations can be found even for large system 
sizes and strong disorder. We illustrated this behavior by considering the sit- 
uation of a two-step global failure. In this particular case, several quantities 
have a discrete jump at the critical point, like the branching ratio ( which 
goes from zero to i.e., from zero to one in the limit of A^o ~^ oo- 

The scenario is quite different for the global load sharing FMB. In this case, 
the type of phase transition is not very clear. There are basically two possible 
ways in which such a transition may occur. In a first-order phase transition, 
we find discontinuous behavior in various quantities as we pass through the 
critical point. Contrary to this case, in a second-order (or continuous) phase 
transition, the fluctuations are correlated over all length scales, therefore, not 
only does the correlation length diverges continuously as the critical point is 
approached, but also other quantities show scaling [24]. A particular case is 
the first-order phase transition close to a spinodal-like instability, for which 
scaling relations can be obtained despite that some macroscopic quantities, 
like the elastic modulus, have a discrete, finite jump at the critical point [25]. 
The question then is whether the global load sharing FBM is a case of a 
second-order transition or a first-order spinodal transition. 

The results obtained with the probabilistic approach seem not to support 
the claim of Rcf.[5] that the failure of a fiber bundle model under a global 
load sharing transfer scheme behaves as a first-order phase transition close 
to a spinodal point. There, by simulating models of electric breakdown and 
fracture, the authors presented ample numerical and theoretical evidence of 
several scaling relations and of the discrete jump of the macroscopic proper- 
ties. We have obtained the same scaling relation for the rate of flber failures as 
the critical point is reached, as well as for the avalanche sizes, which also di- 
verge at the transition. Note that the scaling exponent derived from numerical 
simulations of the continuous version of the probabilistic approach fits very 
well the mean field result 7 = |. Besides, the fraction of unbroken fibers just 
before the global rupture has a discontinuity, which is not size dependent for 
big systems as can be observed in Fig. 5. It has also been shown by analytic 
means that the distribution of avalanche sizes in the case of global load shar- 
ing transfer rule follows a universal power law with an exponent —5/2 [26-28]. 
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However, we should notice that in driven disordered systems, the concepts re- 
lated to spinodal micleation are not sufficiently well established. In particular, 
it seems that the finite size of the system makes it possible that the last catas- 
trophic avalanche provokes the elastic modulus to drop to zero with a finite, 
and very likely, N-dependent jump. Thus, from our point of view, that is not 
enough to set the conclusion that fracture can be described as a first-order 
phase transition. 

Our alternative point of view is to consider the branching ratio defined pre- 
viously as a measure of the distance of the system from the critical point. 
According to the results obtained, the branching ratio goes continuously from 
zero to one. The branching ratio measures the probability that a fiber that 
has failed triggers none, one or more breaking events. Therefore, a branching 
ratio equal to one at the critical point implies that, on average, one failure 
will induce another failure so that the system reaches a state where any per- 
turbation propagates across the entire system, which is, in essence, a critical 
phase. This invokes a continuous phase transition as claimed in other analysis 
of fracture models [29]. Note, additionally, that what changes discontinuously 
at (7c is the rate of change of ( rather than ( itself. 

In summary, we have introduced a probabilistic approach to fiber bundle mod- 
els which allows to avoid the fiuctuations near the critical point. The different 
rupture behaviors of the system for the global load sharing scheme and the 
local one were shown. For the short-range interactions case, the rupture is 
no-doubt of the first-order phase transition type. For the global transfer rule, 
several scaling relations were obtained with mean field critical exponents. The 
branching ratio was defined as an appropriate order parameter. According to 
the results obtained, the branching ratio goes continuously from zero to one. 
This suggests that fracture in heterogeneous systems with long range interac- 
tions can be described as a phase transition of the second-order type, at least 
within the FBM picture. 



5 ACKNOWLEDGMENTS 

Y.M thanks A. Vespignani and H. J. Jensen for very useful and stimulating 
discussions. Y. M also thanks The Abdus Salam International Centre for The- 
oretical Physics for hospitality and financial support. This work was supported 
in part by the Spanish DGICYT under Project PB98-1594. 



References 



13 



[I] A. Garciamartin, A. Guarino, L. Bellon and S. Ciliberto, Phys. Rev. Lett. 79, 
3202 (1997). 

[2] A. Petri, G. Paparo, a. Vespignani, A. Alippi and M. Constantini, Phys. Rev. 
Lett. 73, 3423 (1994). 

[3] Statistical Physics of Fracture and Breakdown in Disordered Systems. B. 
K. Chakrabarti and L. G. Benguigui, Clarendon Press, Oxford (1997), and 
references therein. 

[4] Statistical Models for the Fracture of Disordered Media. Editors, H.J. Herrman 
and S. Roux, North Holland (1990); M. Sahimi, Phys. Rep. 306 (1998) 213, 
and references therein. 

[5] S. Zapperi, P. Ray, H. E. Stanley and A. Vespignani, Phys. Rev. Lett. 78(1997) 
1408; Phys. Rev. E. 59, 5049 (1999). 

[6] J. V. Andersen, D. Sornette, and K.-T. Leung, Phys. Rev. Lett. 78, 2140 

(1997) . 

[7] B.D. Coleman, J. Appl. Phys. 29, 968 (1958) 
[8] H.E. Daniels, Proc. Roy. Soc. A183, 404 (1945). 

[9] Y. Moreno, J. B. Gomez, A. F. Pacheco, Phys. Rev. Lett. 85, 2865 (2000). 

[10] M. Vazquez-Prada, J. B. Gomez, Y. Moreno and A. F. Pacheco, Phys. Rev. 
E. 60, (1999) and references therein. 

[II] D.G. Harlow and S.L. Phoenix, J. Composite Mater. 12, 195 (1978). 
[12] F. Kun, S. Zapperi, and H. J. Herrmann, preprint |cond-mat / 9908226] . 
[13] Y. Moreno, J. B. Gomez, and A. F. Pacheco, Physica A 274, 400 (1999). 
[14] R. da Silveira, Phys. Rev. Lett. 80, 3157 (1998). 

[15] D. Sornette, J. Phys. A 22, L243 (1989); J. Galambos, The Asymptotic Theory 
of Extreme Order Statistics (New York, Wiley, 1978). 

[16] A. Hansen, P. C. Hemmer, Physics Letters A 184, 394 (1994). 

[17] V. I. Raisanen, M. J. Alava, and R. M. Nieminen, Phys. Rev. B 58, 14288 

(1998) . 

[18] H. J. Jensen, Self- Organized Criticality (Cambridge University Press, 1998), 
and references therein. 

[19] G. Caldarelh, C. Castellanos, A. Petri, Physica A 270, 15 (1999). 

[20] T. E. Harris, The Theory of Branching Processes (Berlin: Springer- Verlag, 
1963). 

[21] O. Sotolongo-Costa, Y. Moreno, J. J. Llovera, J. C. Antoranz, Phys. Rev. Lett. 
76, 42 (1996). 



14 



[22] J. B. Gomez, D. Ifiiguez, and A. F. Pacheco, Phys. Rev. Lett. 71, 380 (1993). 

[23] B. Q. Wu, P. L. Leath, Phys. Rev. B. 59, 4002 (1999). 

[24] J. Cardy, Scaling and Renormalization in Statistical Physics, Cambridge 
Lectures Notes in Physics 5, Cambridge (1996). 

[25] L. Monette, Int. J. Mod. Phys. B 8, 1417 (1994). 

[26] P.C. Hemmer and A. Hansen, J. Appl. Mech. 59, 909 (1992). 

[27] M. Kloster, A. Hansen, and P. C. Hemmer, Phys. Rev. E 56, 2615 (1997). 

[28] D. Sornette, J. Phys. A 22, L243 (1989); J. Phys. I (Prance) 2, 2089 (1992). 

[29] D. Sornette, J. V. Andersen, Eur. Phys. J. B, 1, 353 (1998). 



15 



Fig. 1. Rate of fiber failure as a function of the distance to the critical point for 
the continuous version of the probabilistic global load sharing model for a system 
of No = 50000 elements and p = 2. A line with the mean-field value 7 = ^ of the 
exponent is plotted for reference. 



Fig. 2. Evolution of the branching ratio as the critical point is approached in the 
continuous probabilistic method {Nq = 50000). Note that at the critical point the 
branching ratio reaches the unity. 



Fig. 3. Behavior of the branching ratio near the critical point. The straight line 
satisfies the scafing relation (1 — C) ~ {ac — cr)^, with /3 = 0.5. This plot has been 
obtained for a uniform distribution of fiber's strengths. 



Fig. 4. Branching ratio for different system sizes (p = 2) for the global load sharing 
case. As the size of the system is increased the value of the branching ratio for the 
last non-catastrophic avalanche approaches unity. 



Fig. 5. Fraction of unbroken fibers just before the final breakdown in the probabilistic 
model of the global load sharing case as a function of the system size (p = 2). The 
horizontal line shows the exact value in the thermodynamic limit. 
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